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Variations in distinct restricted spaces of wave functions generate distinct density functionals 
(DF)'s. In particular, angular momentum projected Slater determinants define a new DF, compat- 
ible simultaneously with angular momentum zero and a mean field description. 
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Consider a finite number of identical particles, fermions for instance. Let a], and a r be their creation and annihilation 
operators at position r. Their number operator reads, A = J dr a\a r . Consider now a nuclear or atomic Hamiltonian, 
H = T + V + U, where A, T = J^ti U, V = £<> 3 =i % and U = £~ 

^ Ui are the particle number, kinetic energy, 
two-body interaction energy and one-body potential energy, respectively. It is understood in the following that u is 
local, (r\u\r') — u(r) 8{r — r'), although generalizations of the theory to non local potentials are available It is 
also understood that u is a scalar, like t and v, although non scalar it's will be reinstated at the end of this note. 

There is some frustration expressed in the litterature 0] at a contradiction between the existence theorem provided 
by Hohenberg and Kohn (HK) which deals with exact ground states having good quantum numbers, and the results 
given by those functionals used in practice, which often provide solutions showing spontaneously broken symmetries. 
' The difficulty with such practical functionals comes, obviously, from their close connection with mean field theories, 
such as, for instance the Hartee-Fock (HF) variational principle. Moreover, the Kohn-Sham (KS) theorem jjj, uses, 
orthogonal orbitals of independent particles, although its driving mean field does include, in principle, all the influence 
of correlations in a complicated wave function. 

That mean field approximations lead to broken symmetries is actually a very good feature; there is no doubt indeed, 
£f~) | for instance, that Ne 20 is a strongly deformed nucleus! But the subsequent projection of a good angular momentum 
■ is mandatory; indeed the ground state of Ne 20 is a + and the relevant density is spherical, not prolate. 

This note proves an existence theorem for a functional that allows simultaneously for mean fields breaking symme- 
tries and for densities reflecting symmetry restoration. For the sake of generality, we can follow Mermin [§J and use a 
finite temperature formalism. In Fock space, with a grand canonical ensemble, this means using the second quantized 
forms of A and H, naturally. We first consider the domain T> of all density matrices D in Fock space, under the 
obvious normalization constraint of a unit trace. Any approximation such as Hartree-Bogoliubov, BCS, shell model 
mixtures, generator coordinate mixtures, finite temperature HF, etc. then makes a restriction 1Z of the domain T>. 
The usual variational principle, with a temperature T and a chemical potential /i, reads, 

£ ■ TrD (H — uA + T InD) 

F M = Min De2 , F, F = i -f- >- . 1 

1 r D 

With unrestricted density operators, it selects the well known, unique position, Dm = exp [— (H — p A)/T] /Z, with 
Z = Tr exp [— (H — p A)/T] , where a smooth minimum, Fm, is reached. Any variation <5D away from Dm creates a 
strictly positive increase of F at second order with respect to ST). The functional double derivative, S 2 F/(ST) <5D'), is 
strictly positive definite at that position Dm in T>. Accordingly, in any, whether local or global, system of coordinates 
for T>, the matrix representing 5 2 F/(5'D ST)') will be invertible. 

From the very definition of Dm, the first order functional derivative, SF/ST), vanishes at that position. Hence, if 
one introduces an infinitesimal variation Su, triggering an infinitesimal displacement of Dm, the only contribution to 
6Fm comes from 8u. In short, SFm = TV Dm <5U. There is no contribution from <5Dm- 

Define the one-body density matrix in coordinate representation, v(r,r') = Tr Dm <A- a r' ■ Its diagonal, p(r) — 
v(r,r), is the usual density deduced from Dm by integrating out all particles but one. Hence J dr p(r) = TrDif A. 
Note also the result, 5Fm — J dr p(r) 8u{r). In other words, 5Fm/5u = p. 

Freeze t and v and consider Fm as a functional of u alone. The Hohenberg-Kohn-Mermin (HKM) process then 
consists in a Legendre tranform of Fm, based upon the essential result, 

This Legendre transform involves two steps, namely i) subtract from Fm the functional product of u and 8Fm/5u, 
namely the integral, J dru(r) p(r), leaving 

Tm = Tr Dm (T + V- ftA + T In Dm) , (3) 
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then ii) consider p as the primary "variable" rather than u, namely consider Tm as a functional of p. Step ii) is made 
possible by the one-to-one correspondence between u and p, under precautions such as the exclusion of those trivial 
variations Su which modify u by a constant only. (For some discussion of more precautions, see for instance [(J and 
0.) The one-to-one correspondence is proven by an argument ab absurdo : if distinct potentials u and u' generated 
(distinct!) Dm and D' M with the same p, then two contradictory, strict inequalities would occur, namely, 

dr [u(r) - u'(r) ] p(r) < F M - F M and / dr [u(r) - u'(r) ] p(r) > F M - F M . (4) 



(5) 



An inverse Legendre transform, returning from Tm to Fm, is then available, according to the property, 

dp 



Now, the same properties and process hold if we assume that, when D is restricted to a smaller domain 1Z, there is 
still a unique position D m for a smooth again, strict minimum, F m = Minoe-R. F. The subscript M should be replaced 
by the subscript m, but otherwise, it is trivial to recover Eq. (J2J and the contradiction Q to implement a one-to-one 
correspondence and a Legendre transform. 

This will be illustrated by the special case of angular momentum restoration, at zero temperature. Limits at T = 
are easy and do not need much discussion. Consider therefore a canonical situation, where is no need for p\ there is a 
well defined particle number A, and density operators are dyadics, D = |?/>)(7/>|, where \ip) is just an A-particle wave 
function. The variational principle, Eq. (1), trivially reduces to the Rayleigh-Ritz form, 

F M = Min,F, f-^, (6) 

and the minimum with unrestricted wave functions generates the exact ground state with its exact eigenvalue Fm, 
parametrized by u. Restrict now the wave functions to be angular momentum projected Slater determinants, 

* = J dflR(n) (7) 

where fl and R label rotation angles and rotation operators, respectively, and \</>) is an unrestricted, hence most 
often deformed, Slater determinant. Angular momentum has been chosen here because every even-even nucleus has 
such a ground state. In a short, transparent notation, |^ ) = Po |</>); where P can be normalized to be idempotent, 
Pq = Po. It commutes with H, a scalar, and we use these commutation and idempotence in the following. 
The key question is now whether there is a unique fo™ for which the minimum, 

(*E |H|tt ) 

F ro -Mm* {Wo) , (8) 

is reached. Notice that, in any case, the restriction of trial functions enforces an inequality, F m > Fm, and that this 
inequality is generally strict, F m > Fm ; it demands very special Hamiltonians and potentials for an equality to occur. 
For deformed nuclei the same minimum F m , which also reads, 

MH P O |0) 

Fm-Mi^ , (9) 

occurs for many choices of cj) . All such Slater determinants are deduced from one another by rotations. Hence they 
generate the same projected state ■ As illustrated by Figure 1, that breach of uniqueness, namely the rotational 
degeneracy, caused at the <f> level by deformation, is completely corrected at the level. There might be further, 
accidental causes of degeneracies of the minimum in the space, but it is clear that the broken symmetry degeneracy 
has been removed by the projector Po . 

It is thus reasonable to assume that, provided all those symmetries broken by HF and similar approximations have 
been restored, a search for a minimum energy, after the restoration, will suffer from no degeneracy any more. 

Under this assumption of uniqueness of the minimum, consider the space of those V^o that are square normalized 
to unity. (This technicality spares us denominators.) Let $ora and be those states which provide F m and F^ 
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Slater determinants 



FIG. 1: One-to-one correspondence u «-» p. A given potential u (right) induces degenerate Slater determinant variational 
solutions (center) resulting in the same density p after angular momentum projection (left). 



for H and H' driven by distinct u and u', respectively. Then, the same argument ab absurdo as that for Q, with the 
same precaution of a non trivial difference between u and u', 



J dr [u(r) - u'{r)\ p(r) = <* 0m | (H - H') |* 0m ) < F m - , 

J dr [«(r) - u'(r)] p(r) - (*^| (H - H') |^ m ) >F m -F^. (10) 

proves that fom and fp m cannot induce the same p. Hence the one-to-one map, u <-> p, holds for the restricted 
variation space. It must be stressed here that the density p is that of the symmetry projected state. It can differ from 
the density of every parent state <fi through the correlations brought by the coherent summation over angles, Eq. JJJ. 

The "smoothness" condition of F around its minimum is necessary to allow for functional derivatives 6F m /5u, 
5 2 F m /(Su5u'), etc. Since there are are no singularities in a space of Slater determinants (although this space is 
curved), and since that subset of determinants which might create singularities because of a vanishing matrix clement 
(</>|Po|</>) has obviously a zero measure, we recover all the properties of the unrestricted variational space. In particular, 
5F m /5u = p. A Legendre transform is possible. Hence there exists a functional of a spherical p, with value 

^m = (*0 m |(T + V)|* m). (11) 

It may be stressed here how a spherically invariant u is in a one-to-one correspondence with a spherically invariant p. 
This is reinsuring for the consistency of this existence theorem, deduced from a restricted variational space. 

Since $|j ra can only be an approximate ground state, then T m > Tm-, compare Eq. J3J without p at T = and 
Eq. Hll|) . We are thus dealing with a new functional. It can only give upper bounds to bound state energies, but this 
is expected for restricted variational spaces, obviously. 

For the case of nuclear forces, with their very repulsive core and need for short range correlations, nothing prevents 
us, formally at least, from refining the variational space by introducing a "correlator operator" a la Jastrow for 
instance, that crushes the wave function at short relative distances between particles, 

F „ =Mi „, <fggZf , _ Jjl*(|«-r,|), 
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where x vanishes at the origin and becomes 1 at large distances. Since C is a scalar, it commutes with Pq. 

Our considerations can be easily generalized to finite temperatures, and to other symmetry restorations than angular 
momentum projections, and to other restricted spaces. It can be concluded that every restricted space and every 
symmetry restoration creates its own density functional, provided the conditions of uniqueness and smoothness of the 
variational minimum are fulfilled. For the discussion of differentiability and fine topological properties of the u- and p 
spaces we refer again to 6]. Up to our understanding, the validity of our existence theorems for symmetry conserving 
functionals is not compromised. 

Such existence theorems, though, suffer from the usual plague of the field: constructive algorithms are missing. 

Another question must now be raised: can one define a KS theory compatible with symmetry restoration? Consider 
again the case of angular momentum 0, for non degenerate ground states of scalar Hamiltonians, with p scalar, as are 
t, v and u. One may then claim the exact Fm and its partner Tm are also scalar. Accordingly, all the tools of the KS 
theory, including the "under-rug-sweeping" exchange-correlation potential v xc (r), must be rotationally invariant. For 
the sake of simplicity, let us forget spin-orbit complications. Then the KS equations read in coordinate space only, 



t + 



dr' v(r — r') p(r') + v xc (r) + u(r) 



<ptm(r) = 0, 



(13) 



Can they be assumed to generate the same density as that of the correlated ground state? 

In Eqs. H 13p . the magnetic degeneracy of the orbital eigenvalues £i is explicit. A complete filling of a core, made of 
all lowest shells but one, and a homogeneous, partial filling of a cloud for the residual particle number z in the next, 
partly filled shell, seem to be in order. It is trivial indeed that the following density, 



p(r) 



E 

i^core 



2L + 1 



E 

MS cloud 



(14) 



is rotationally invariant. Here L is the orbital label for the cloud, that first shell above the filled core, and M labels 
the corresponding degenerate orbitals. Can such a density be related to the angular momentum projection of any 
state in the subspace spanned by all the degenerate Slater determinants available when filling z among the degenerate 
2L + 1 orbitals? Notice, incidentally, that necessarily z ^ 1 and z ^ 2L, since such situations would only allow total 
angular momentum L. Notice also that particle- hole symmetry that relates the coupling of angular momenta in a 
situation with 2L+1-Z to that with z. For simplicity, we thus set 2 < z < L. Then we forget about the core and 
investigate the densities of those angular momentum projected determinants made of mixed orbitals inside the cloud. 
By mixed orbitals, we mean arbitrary linear superpositions of the magnetic numbers M; such mixtures play the role 
of deformations. What is the density of such projected "deformed" determinants? 

Consider L — 6, z — 4 for an illustrative example. Standard angular momentum techniques prove that one can then 
make two distinct states, ^>oa and with total angular momentum 0. Therefore the projected Slater determinants 
must be linear combinations a ^>o a + (3 ^ob, with densities a 2 p a + ft 2 Pb + 2 a (3 p a b, where p a b is the "transition density" , 



Pab{r) = (*ob| 



at a r 



*0a>, 



(15) 



while p a and pb are the densities of \& a and ^h, respectively. The latter densities are the same, however. Indeed, 
since all orbitals have the same radial shape, <p (|r|) , there are only angular integrals to consider when calculating the 
densities. For instance, 



Pa(r) = 4</?(|r|) 2 / df 2 df 3 df 4 \ip a (f, f 2 ,f 3 ,f 4 ) | z 



(16) 



where fj designates the polar angles of a particle , and ip a carries the angular properties of ^oa- Because of the scalar 
nature of p a , it is obvious that the angular integral gives a result, p a (r) = 4ip(\r\) 2 , which does not depend on the 
singled out f. The same holds for ^ob- Moreover, the transition density reads, 



pab(r) = 4^(|r|) 2 df 2 df 3 df4 ^6 (r, f 2 , r 3 , f 4 )* ipa{f 1 f 2l f 3n f i ) 



(17) 



and this must again reduce to a constant with respect to f because of the rotation invariance of both \l/oa and ^/q . 
The constant, however, must vanish. This is because, upon an integration with respect to ?*, one must recover the 
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orthogonality of ^Q a and ^06- It can be concluded that every "deformed" determinant made of orbitals taken from 
the spherical "cloud" shell has the same density. In Eq. I|14|). the sometimes called "filling" approximation term, 



2L + 1 



<p(\r\) 



(18) 



M dcloud 



is not an approximation; it actually makes the exact ansatz for an exact, projected density. 

The burden of an accurate description of p then falls upon the radial shape y>(|r|), that is the output of the KS 
equations, Eqs. (|13|l . This makes it all the more important to better describe the often badly known potential v xc . 

Generalizations to other values of L and z and to more than two independent states with total angular momentum 
zero are obvious. All determinants from the cloud shell give the same density after projection. 

To conclude, a word is in order about non spherical potentials u, and more generally symmetry breaking potentials 
u. As long as one is interested in Eq. ©, or rather in that proper form that takes into account the fact that u does 
not commute any more with the symmetry restoring projector P, (still normalized to be idempotent,) 



Mm, 



IPHP 



(19) 



it is clear that only that symmetric part, PUP, is activated. Hence a non symmetric u brings nothing. The situation 
is different for the KS theory. Indeed, both the Hartree-like field, / dr' v(r — r') p(r'), and a symmetry violating u lift, 
in general, the unlikely degeneracies of the KS orbitals, and an unprojected density results from, 



t + / dr' v(r — r') p(r') + v xc (r) + u{r) 



<Pi(r)=0, p(r) = Y \<Pi{r)\ 



(20) 



The sum, J2i > 1S taken, naturally, upon the lowest orbitals, except in very few anomalous cases. Consider a limit 
where u becomes symmetric. As is well known, there is no reason why, at that limit, a symmetric state should emerge 
out of Eqs. H20[l ■ One needs an ansatz like Eq. 1|14|) to enforce the symmetry. Under such an ansatz, which should 
be tuned for the partial filling of quasi degenerate orbitals, as with a finite temperature, the Hartree potential and 
a slightly symmetry violating u will indeed create a slight degeneracy lifting only. There is thus a serious difference 
between an unconstrained KS approach and the symmetry constrained one. 

To summarize this note, two distinct results are claimed, namely i) existence theorems for density functionals 
explicitly enforcing a symmetry by means of projectors chiseling variational spaces, and ii) a symmetry respectful 
version of the KS theory. The latter comes from an adjustment of the definition of the density; KS orbitals need to 
be incoherently mixed to account for symmetry projection. 
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